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A class of the symmetrically frustrated SU(iV) models is constructed for quantum magnets based 
on the generators of SU(JV) group. The total Hamiltonian lacks SU(iV) symmtry. A mean field 
theory in the quasi-particle representation is developed for spin liquid states. Numerical solutions in 
two dimension indicate that the ground states are gapless and the quasi-particles are Dirac particles. 
The mechanism may be helpful in exploring the spin liquid phases in the spin-1 bilinear-biquadratic 
model and the spin-orbital model in higher dimensions. 

PACS numbers: 75.10.-b 

I. INTRODUCTION 

Exploring spin liquid has been a recurrent theme in the research of strongly correlated electron systems over the 
past decades [HEHSEIIHII. The picturc for s t 

>in liquid is based on the concept of resonating valence bonds 
(RVB), which says all spins form spin singlet pairs [sL floL fllj - Spin singlet pairs consist of either nearest neighboring 
spins or long-range separated spins. Quantum coherence of spin singlets determines the quantum properties of the 
state. It is realized that a short-range RVB state usually exhibits a finite gap for spin excitations [lj], while a long- 
range RVB state possesses antiferromagnetic long-range order (LRO) [l3l[l3|. Nowadays a lot of systems are found 
to exhibit the behaviors of spin liquid. One route to explore spin liquid is to focus on geometrically frustrated spin 
systems or systems with breaking translational invariance. The formers include the Kagome lattice model and the 
Ji — J2 model. The latters include the spin Peierls systems and the plaquette RVB systems. In these systems the 
quantum and/or geometric frustration is the key point to realize spin liquid. Another route is to deal with a series of 
low-dimensional spin systems with Haldane gap or the spin systems with orbital degeneracy which enhances quantum 
frustration. Such systems neither break the translational invariance nor exhibit geometrical frustration. 

In this paper we propose a class of quantum frustrated SU(iV) models for quantum magnets, and develop a mean 
field theory for the spin liquid involved in it. The models are constructed by the SU(iV) generators, but has no 
SU(iV) symmetry because of the competition between the generators in two different representations. For small N 
it is shown that these models are equivalent to several physical models which have been extensively studied, and fall 
into the same mathematical structures in the particle representation. To explore the spin liquid a mean field theory 
based on the RVB state is developed. In two dimension (2D) it is found that the ground states for N = 3 and 4 
are gapless spin liquid and the collective excitations are Dirac particles. As examples we apply the theory to spin 1 
bilinear-biquadratic system and spin 1/2 systems with double orbital degeneracy. 

II. MODELS 

Let's start from a class of Hamiltonian 

h = Ji £ J?(n)j;(r 3 )-J 2 J2 WnWfa), (i) 

which is defined on a lattice. J/f(r 2 ) are the N 2 — 1 generators of SU(iV) group, and satisfy the algebra 

[Jpin), J£(r,-)] = S tJ (#J£(ri) - Sp:(r t )) . (2) 

We choose Jjffa) as the fundamental representation with a single box in Young tableau. J\ and J2 are two coupling 
constants. 

The first term in Eq.JTJ possesses the SU(AT) symmetry because the operator, Jv( r i)J}l( r j) = Aj> serves as the 
permutation operator, which swaps two quantum states at sites i and j. This term in various forms has been studied 
extensively. The case for N = 2 is the well known Heisenberg model, which was solved exactly and exhibits a gapless 
ground state in one dimension (II?) ^5|. While in higher dimensions it is well established that the ground state 
possesses antiferromagnetic LRO on hypercubic lattices. The case for N = 3 was known as the exchange model for 
ferromagnetism ^| 0, . For N = 4 it was proposed that a spin liquid exists on the triangular lattice |2(], El , 
but a Schwinger boson theory gives a state with LRO j2^| ■ 
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The second term in Eq.Q) also has the SXJ(N) symmetry on a bip artite lattice, which means that J2 connects two 
lattice sites belonging to two different sublattices [23|, |2J, |2{| |2(| |27], |2g . In this case the generators at two sublattice 
sites are expressed in conjugate representations since J£(rj) = [J£ (rj)] , and thus the SU(iV) symmetry survives. In 
ID the ground state is gaped for N = 4. In 2D, the ground states are found to exhibit the Neel-type LRO with 
broken SU(JV) symmetry for N < 4, and possibly to be a spin liquid for N > 4 I2I I27L |28| . 

The mixture of the two terms makes the whole Hamiltonian in Eq. JJ) deviate from the two SU(iV) symmetries. 
We propose that this symmetric frustration may also lead to spin liquid. To study the phase diagram of the model in 
Eq. Q we introduce a set of creation and annihilation operators to rewrite the Hamiltonian in the second quantization 
representation. In the SU(JV) representation each site has N quantum states \fi) so that we may introduce N pairs of 
operators &| and b itl : \i,fi) — b\ |0) with the vacuum state |0) at site i. In this way we can construct the operator, 

J„(ri) = bl v bin, with a constraint for single occupancy, X)u=i = 1j 011 eacn site. This is the so-called hard-core 

condition even if the particles are bosons. Interestingly it is found that the generators satisfy the SU(iV) algebra in 
Eq. (2) for either Boson or Fermion representation. So the Hamiltonian is reduced to 

h = j iJ2 p * J 2 E B h B v + E x 4 E b l b ^ 1 ) ■ (3) 

ij ij i \ IX / 

where the bond pairing operator Bij = bj^bi^ and the Lagrangian multipliers Xi are introduced to realize the 
constraint of single occupancy. The permutation operator can be expressed as Pij = X^i/ bl^bi^bj^bj^ = — q + qFj-Fij 
with Fij = "jpbifj,, and q = 1 for bosons and —1 for fermions. Several physical systems are shown to belong to 
this Hamiltonian. It has been known for a long time that the spin SU(2) operators with S can be used to express 
the SU(iV) generators. Schrodinger developed an expression for the permutation operator in terms of (Sj • Sj) m for 
m = 0, 1, ■ • ■ , 2S |29l ]. (Please see in Appendix A.) We present several examples which have the equivalent form in 
cither the Fermion or Boson representation. 

For N = 2 the model in Eq. is equivalent to the spin 1/2 XX Z model, 

// E [J**SfSf + J**S?S* + J ZZ S*S]} . (4) 

ij 

By defining b\ |0) and b\ |0) as the two eigenstates of the operator Sf , the spin operators Sf can be expressed in terms 
of b operators, 

(5a) 
(5b) 
(5c) 

In this way, one can show that the model has the form of Eq. for N = 2 with Ji = (J zz + J xx )/A and J 2 = 
(J zz — J X k)/4. The model with a similar form was studied by Leone and Zimanyi |3(j . When J zz = J xx , it returns to 
the well-known Hciscnberg model. 

For N = 3 it is the spin 1 bilinear-biquadratic model, 

H = [cos0 S, • Sj + sm<jy (Si ■ S^) 2 ] , (6) 

ij 

which is one of the prototype models exhibiting Haldane gap in ID antiferromagnet |3ll l32|. The phase diagram 
for ID is well established that the energy gap persists in the wide range, — tt/4 < <j> < n/4 |33l 13^. l35l |. The phase 
diagram for 2D and 3D was studied by means of quantum Monte-Carlo p3q . Non-zero quadrupole moment at zero 
temperatures is found at the region, — n < <p < 0. For spin 1, each site has three states \nii) with m, = — 1,0, +1 
(i = 1, 2, 3) according to the eigenvalues of Sf . We reorganize the three states and define three operators, 

&1|0) = -Ldm^ + K)), (7a) 
bt\0) - \m 2 ), (7b) 
6||0) - -L(K)-|m 3 )). (7c) 
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In terms of b operators, the three spin operators can be written as 

SI 
Sf 

Sf 



f = i(bl 2 b a 


-bib 




-b\ 2 b 




-b\ 3 b 



(8a) 
(8b) 
(8c) 



On this basis the Hamiltonian Eq. JBJl is reduced to the form of Eq. © for N = 3 with J\ = cos0 and J2 = 
ypl cos ((f) + 7r/4) 0. At the points, (j) = — 37r/4, 7r/4, ±7r/2 the model possesses the SU(3) symmetry 15]. It is 
solvable at tan0 = 1/3 and the ground state is a valence bond solid [3?i| . 

For TV = 4, we can present a model in terms of operators S = 3/2 just like the cases for N = 2 and 3. Alternatively 
we present the spin-orbital model [38j . 



H = (Si • S 3 - + T. t • T,) + 4V J2 & ■ ( T * ' T j) 



(9) 



with spin S — 1/2 and orbital T = 1/2, where both of the operators, S and T, satisfy the SU(2) algebra. This model 
was used extensively to describe the orbital physics in transition metal oxides |39j| . At J — V the model has the 
SU(4) symmetry. There are four simultaneous eigcnstates for spin S z and T z at each site, \Sf = ±|,Tf = ±±). So 
we introduce four creation operators as follows, 



IK* -£( 

4|o> = 



6 10) = 



1 

3 



1 1 

2'~2 
1 1 
2'~2 
1 1 

-2' + 2 
1 1 

~2' + 2 



1 1 

2' + 2 
1 1 
2' + 2 
1 1 
2'~2 
1 1 
2'~2 



(10a) 
(10b) 
(10c) 
(lOd) 



By introducing a four component spinor ^ = \ b\, b\, b\j , we can express the spin and orbital operators in terms 
of * t and 



s x = 


i* f (<r y ® cr ) *, 


(11a) 


s v = 


(<r x ® o-j,) *, 


(lib) 


s z = 


i* t (-er,®o- tf )*, 


(11c) 


rpX 


i* f (-er w ® (7,) *, 


(lid) 


T y = 


i* f (-<r w (8 (T X ) *, 


(lie) 




(cr <g> cr y ) 


(llf) 



Once again one can show that the spin-orbital model falls into the general form of Hamiltonian as Eq. @ with 
Ji = (J+ V)/2 and J 2 = (J - V)/2. 

So far we have shown that the three physical models can be expressed in a general form in terms of SU(iV) generators 
as Eq.(l) or in the particle representation as Eq. (2). For larger N, the models are related to some high spin systems 
with orbital degeneracy, for example, the spin-1 system with double orbital degeneracy |4*cj . 



III. MEAN FIELD THEORY 



Now we turn to the phase diagram of the model on a d-dimensional hypercubic lattice. We may choose the operators, 
6's, as either bosons, [6^, &tj = 8^8^, or fcrmions, j&j^, b^ v X = 5^8^. In principle the bosons tend to condensate to 
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the lowest energy state at low temperatures and to form a quantum ordered state, while the fcrmions tend to form a 
Fermi sea and a quantum disordered states. In practice it is inevitable to introduce adequate approximation schemes 
to deal with a many-body system in either the Bose or Fcrmion representation. Since our purpose is to search the 
regime of the spin liquid phase, we choose the Fermion representation in this paper. Without loss of generality we 
can assume that the system is defined on a rf-dimensional simple cubic lattice, and take J±, J 2 > 0. In this way the 
Hamiltonian is semi-negative in the Fermi representation. Now the second term in Eq. Q plays a role of an attractive 
interaction between pairing valence bonds. [2q The Hubbard-Stratonovich transformation is performed to decouple 
the Hamiltonian into a bilinear form. Two types of mean fields arc introduced, 



Ai(fc) = JiJ2( F ^+s)e lk ' S = 2 J ± FJ2 cos k a ; 



A 2 (fc) = J 2 J2(Bi,i+ 



s) e 



ik-5 



2J 2 B sin k a 



(12a) 
(12b) 



where F = (Fij+s) = (Fij-s) and B = i (Bij+s) = —i(Bi,i—s)i and ol = x, y, ■ ■ ■ . The bracket (•••) represents 
the thermodynamic average. The chemical potential Aj is taken to be site- independent, \ = A, which can be also 
regarded as a mean field. In the momentum space the mean field Hamiltonian is 



b 1 



fe , [a k 

—XNa + dN A JiF 2 + dN A J 2 B 2 . 



(13) 



where e(k) = A — Ai(fc), Na is the total number of lattice sites. This is a Bardccn-Cooper-Schrieffer (BCS) type of 
Hamiltonian for iV-component degenerate fcrmions. By performing the Bogoliubov transformation, 



7fcp = Ukbkf, - v k bl kll ] jl kfJ = u k bl kfi + v k b kll 



with the coherence factors satisfying 



1 + 
1 - 



e(k) 



Lj{k) 



2u k v k 



A 2 (fc) 



one can diagonalize the Hamiltonian as 



where the spectrum and the ground energy are 



(14) 

(15a) 
(15b) 
(15c) 

(16) 
(17) 



E = ~J2"( k ) + ^y^AA^A + dN A ,hF 2 + dN A J 2 B 2 . 



(18) 



The spectra of ui(k) are A^-fold degenerate for the quasiparticles. From the free energy for the fermion gas 

fi = -^E ln (l + e-^ (fc) )+i?o, 
P , 



(19) 
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we obtain a set of the mean field equations by optimizing the free energy with respect to the mean fields F, B, and 
A, 

dk e(k) Buj(k) N-2 , N 

dk e W (-E„cosM tanh M/0 = MF 



(27r) d w(fc) 2 TV 

dfc A 2 (fc)(E a smfcq) tanh M fc ) = 2dB (2Qc) 



(2^) d 2 AT 

Thus the mean field Hamiltonian is solved together with the self-consistent equations for the three types of mean 
fields. 

We have N degenerate spectra for quasi-fcrmions. These arc not the collective modes, and cannot be measured 
explicitly. In order to calculate spin susceptibility we define the Matsubara Green's function in the form of a 2 x 2 
matrix, 

rrir^- It ( M^L(O) Mr)&_ fc „(0) U , 9U 



where the bracket (■ • ■ ) means that the thermodynamic average is made, and the factor T T is the imaginary time 
order operator. The Green's function of the frequency is obtained as follows, 

n ru 4 \ iuj n a + e(k)a z + A 2 (k)a x 

G M (fc,iw„) = — -3 — , (22) 

\lU) n ) -UJ z (k) 

where ui n = (2n + 1)tt//3 for all integer n. Due to the degeneracy of the spectra the Green's functions are independent 
of the index /i. The imaginary time dynamic correlation function for Jf} is defined as 

Xfl »(q,iu n )= dTe^ T (T T J»(q,T)Jt;(-q,0)). (23) 



The imaginary part of x^il^) is worked out as 

/dk 
—^Siuj - uj{k) - u;(k + q)) 
{2tt) 

x {u 2 k v 2 k+q + u k v k u k+q v k+q 8^ u ) (24) 

for uj > at T = 0. It is shown that ImXfiu{q,^>) becomes non-zero only when u> > A gap = 2mina;(fc). Thus if the 
spectra for quasiparticles have an energy gap the collective excitation for the dynamic correlation function also have 
a finite energy gap. The energy gap A ffap can be evaluated by solving the mean field equations. 

IV. TWO-DIMENSIONAL GAPLESS SPIN LIQUIDS 

It is easy to see that when J2 = or B = the spectra for fermions are reduced to an ideal fermion gas. In the 
case the system is a type of gapless spin liquid as discussed in the SU(2) and SU(4) systems. When J 2 increases a 
non-zero solution for B and F is obtained. We take J\ = Jcos9 and J2 = Jsin9 with < 9 < ir/2. For finite N, 
the mean field equations are solved numerically at T = 0. In one dimensional case a finite gap is found in a large 
regime of the parameter space. In this paper we focus on 2D systems. Numerical solution for non-zero B and F has 
lower energy than other mean field solutions with F — or B = when 8 < 1.075 for N — 3 and < 0.867 for N = 4. 
Opposite to the ID cases non-zero B and F do not produce an gaped phase. Instead they lead to a gapless spin liquid 
with A gap = 0. Notice that along the line = — fc°, we have A 2 (fc) = 0. And on this line if e(k) = the spectrum 
becomes u>{k) = 0. So we get the solution, 

k° x = -k° = ± arccos — — . (25) 
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FIG. 1: Above: the value of k% at two points for the fermi surface (feS,— &2) and (— f° r iV = 3,4 in square lattice. 
Bottom: The velocities of the Dirac quasiparticles along the different directions. For each N there exists a fixed point at which 
the velocity becomes isotropic. The angle 6 = arctan(J2/Ji). 

Numerical values of the parameters show that jj^p < 1 always holds for the mean field solution. \k®\ is plotted in 
Fig. 1. In other words, the Fermi surface actually consists of two points located at fc° = — = ± arccos 4/ ^ F ■ We 
expand the dispersion relation near the Fermi point (k x ,ky) with 

Ak x = k x — = Ak cos ip; (26a) 
Ak y = k y -k° = Ak sin tp, (26b) 

the spectrum becomes linear with respect to Afc, 

w (k° + Ak) ~ v°Ak (27) 

where the anisotropic velocity is 

v° = XV2^Jcl + (c 2 2 - c\) sin 2 (ip + tt/4) (28) 

with c\ = (2JiF/X) 2 — 1/4 andc?, = B 2 tan 2 9/F 2 . The dispersion relation is linear and the quasiparticles are type of 
Dirac particles. Its velocity is in general anisotropic and has its minimal for ip = 7r/4 and its maximal for — 7r/4 for a 
specific 9. At 9 = 0.951 for N = 3 and 0.786 for N = 4 the velocity becomes isotropic. Beyond the regime in Fig. 1 
we do not anticipate the present mean field theory is still valid because it is already known that the long-range order 
exists in the ground state in 2D when J\ = H^, |2j| . 

We have also calculated the static correlation function Xnn an d Xnv (/•* 7^ v) for N = 3 and 4. We found that the 
results for the two systems are very similar. So we only present the results for N = 4 here. In Fig. 2, we plot the 
diagonal susceptibility, x^ilxi 9y)> for different couplings. All values of \ are convergent, which indicates the absence 
of long-range order as expected in a spin liquid state. We observe that x has its maximal value at points (±tt, ±tt) for 
small 8 (or Ji), which shows the static antiferromagnetic correlation dominates. The maximal values decrease with 
the increasing 9 (from Fig. 2(a) to 2(b)), and Xw becomes almost flat for large 8. In Fig. 3, we plot the off-diagonal 
Xfiv (A* 7^ v )- We found that the maximal value of x changes from the points (±7r, ±7r) to (0,0) with 9 increasing. 
Thus the antiferromagnetic correlation decreases very quickly, and the ferromagnetic correlation dominates at large 
9. For an intermediate 9 either the diagonal or the off-diagonal x nas a relatively flat structure, which implies that 
the state is a very well-defined spin liquid. 

Though the mean field theory has lots of disadvantages, it is still used extensively to study the many body system 
with strong correlation after some adequate transformations. Until present there has no definite conclusion on whether 
the 2D SU(iV) (N > 2) antiferromagnetic model with J x > and J 2 = possesses LRO ground state [iH 12(1 12I l22j . 
In the present mean field theory it is assumed that the ground state is a spin liquid in the case. Quantum Monte 
Carlo simulation may provide helpful information. The phase diagram of the 2D spin 1 bilinear-biquadratic model 
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FIG. 2: Diagonal susceptibility X* a {qx,q y ) for N = 4 and different couplings 6 = 0.300 (a), 0.500 (b), 0.600 (c), 0.700 (d), 
0.786 (e), 0.850 (f). x-axis: q x 6 (— 7r,7r), y-axis: q y 6 (— tv,tt), z-axis: \aa G (0.00, 0.28). 




FIG. 3: Off-diagonal susceptibility for TV = 4 and different couplings 6 = 0.300 (a), 0.500 (b), 0.600 (c), 0.700 (d), 0.786 
(e), 0.850 (f). The axis labels are the same as in Fig. 2.. 



in Eq. © was studied numerically for — tt < < |36j . It is found that the ground state is antiferromagnetic for 
— 7r/2 < cf> < 0, which corresponds to the regime 7r/4 < 8 < tt/2 for N = 3 in Fig. 1. As there lacks data in the regime 
< 6 < a/A (w 0.785) (0 < cj) < 7r/4) the numerical calculation does not exclude possible existence of spin liquid in 
the case of J\ , Ji > though the antiferromagnetic fluctuations may be very strong. 



V. CONCLUSION 



In conclusion, we have proposed a class of frustrated SXJ(N) models for quantum magnets based on the generators 
of SU(iV) group, and showed that several physical models have the same mathematical structure in the Fermion or 
Boson representation. The model we constructed breaks SU(7V) symmetry when J\ and Ji are non-zero simultaneously, 
although the two parts in the model do have different SU(iV) symmetries, respectively. Usually geometric quantum 
frustration on lattices or coupling structures is one of the main driving forces to generate spin liquid. The mixture of 
two parts with different SU(A r ) symmetries provides us a possibly new way to seek for spin liquid in higher dimensions. 
In this paper a mean field theory with paired valence bonds for quasi-particles is developed to explore this kind of spin 
liquid. A gapless spin liquid exists in a large extent regime in two dimension. It is believed that the mechanism of the 
spin liquid is the competition of two interactions with different SU(A r ) symmetries. Thus the symmetric frustration 
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is a new possible route to explore the spin liquid in higher dimensions. 
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VI. APPENDIX: 

A. Spin exchange operator 

The role of the permutation operator is to swap the two quantum states at two different sites, 

Pij \i,a;j,(3) = \i,(3;j,a) . 

Pfj = 1 and its two eigenvalues are ±1. The connection between the permutation operator and spin operator was first 
established for spin 1/2 by Dirac^J, and then generalized to the case of arbitrary spin 5 by Schrodinger [2S|, who 
showed that 

2S 



p%j — ' A n (Si ■ Sj) 



n=0 

where the coefficients A n typically have the values, for S — 1/2, Aq = 1/2, A\ = 2; for S = 1, Aq = —1, Ai = 1, 
A 2 = 1; for S = 3/2, A = -67/32, A 1 = -9/8, A 2 = 11/18, A 3 = 2/9. The general expression for P v] is 

P tj = [{2S)\r 2 [M-(2S-l)2S} 

x [M - (2S - 2) (25 - 1)] • • • [M - 2] M 
- [(25 - l)!] -2 [M - (25 - 2) (25 - 1)] 
• • • x [M - 2] M 
+ ... + (_l) 2 s- 1 M+(-l) 2S 

2S 



X;(-l) n [( 2 5 - nW 2 [M - (25 - 1 - n) (25 - n)] 



n=0 

• • • x [M - 2] M 



where 



M = (Sj + Sj) 2 . 



B. Projection operators for two spins 



Using Sij = Si ■ Sj = | (Si + Sj) 2 — 5(5 + 1) and the completeness relation for the projection operators 



Ejlo P J («>J) = one nnds 



(s y ) m = E 



2S r 



.7=0 



2 J(J + l)-5(5 + l) 



where Pj(i,j) is defined as the projection operator for total spin J of two spins 5. For spin 1/2, 



For spin 1, 



For spin 3/2, 
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